International Journal of Trend in Scientific Research and Development, Volume 1(4), ISSN: 2456-6470 www.ijtsrd.com 


Optimization of 36cm2 Effective Area on 
Serpentine Flow Channel of PEMFC 

Dr.V. Lakshminarayanan 

Department of Mechanical Engineering, 

B V Raju Institute of Technology, Narsapur, Telangana, India 


ABSTRACT 

The Proton Exchange Membrane Fuel Cell (PEMFC) 
is an electrochemical device and its performance 
depends on the operating and design parameters. In 
this paper, optimization of design parameters like 
flow channel design, number of flow path, channel 
depth and width, cross section of the flow channel and 
operating parameters as operating pressure, 
temperature, relative humidity, mass flow rate of the 
reactant gases and stoichiometric ratio of the reactants 
on serpentine flow channel of 36 cm2effective area of 
the PEMFC was considered. Creo software used for 
modeling and CFD software packages used for 
analysis of PEMFC. The optimization was done bf 1 
Taguchi method using Minitab 17 software. Based on 
the optimization study, the Landing to Channel width 
ratio (L: C) - 1:2 has produced maximum power 
density of 0.276 W/cm2 also square of response factor 
(R2) was achieved by Taguchi method as 96.9%. 
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I. INTRODUCTION 

The environmental impacts of non-renewable power 
sources and energy scarcity in many countries, 
focuses on renewable energy development. The 
Polymer Electrolyte Membrane fuel cells are 
environmentally friendly power source which is 
suitable for powering both portable devices and 
mobile application due to their high energy density 
and lower operating temperature range viz, 30oC - 
70oC.lt is Eco-friendly power source which is suitable 
for poweringboth portable devices and mobile 
application due to theirhigh energy density and lower 
operating temperature range [1]. In this paper the 
performance enhancement of fuel cell by optimizing 
the influence of various operating and design 
parametersusing CFD Fluent 14.0 and MATLAB 
software packages. The PEMFC consists of polymer 


solid electrolyte membrane placed between an anode 
and a cathode. The electrochemical reaction produces 
electric current along with water and heat as 
byproducts. To attain high current, peak power 
density, proper temperature distribution and optimum 
water management the various flow channel design 
was used. S. Simple [2] addressed in his study that the 
performance of PEMFC has been influenced by the 
flow channel path length and the flow field design 
with 1.2,2 stoichiometric ratios, inlet temperatures 
were 70°C and 80°C, operating pressure was l.Olbar 
in the anode and the cathode respectively. 

The numerical analysis on six different cross-sections 
of the flow channel like square, triangle, 
parallelogram 14o, parallelogram 26o, trapezium and 
inverted trapezium having 1.25 cm2 effective area of 
single pass PEMFC were carried out by 
Lakshminarayanan et al [3]. It was concluded that, 
square flow channel had a peak power density of 
1.133 W/cm2. Nicholas S. Siefert [4] concluded his 
study that the serpentine channels are very effective to 
get rid of the liquid because of its high gas velocity. 
The performance enhancement of the combined effect 
of operating and design parameters (operating 
pressure, temperature and inlet mass flow rate of 
reactant gases and rib to channel width ratio) of single 
pass serpentine and interdigitated flow channel with 
25 cm2 active area of PEMFC carried out by 
Lakshminarayanan and Karthikeyan [5]. The results 
revealed that the maximum power density of 
interdigitated flow channel with landing to channel 
width 1:2 showed better performance than the 
serpentine flow channel with same design parameter. 
The various operating parameters like cell 
temperature, pressure, reactants on anode and cathode 
flow rate has been investigated with triangular 
channel geometry on 25 cm2 active area of PEMFC 
by Khazaee et al. [6]. The results showed that an 
increase in the inlet temperature of reactants, cell 
temperature and inlet pressure can enhance cell 
performance of the PEMFC. The effects of 
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interdigitated flow channel with traditional flow 
channel, the effects of the flow area ratio and the 
baffle-blocked position of the interdigitated flow field 
on the performance of PEMFC were examined 
experimentally by Yan et al [7]. The results concluded 
that, the cell performance can be enhanced with an 
increased inlet flow rate of reactant and cathode 
humidification temperature. The increasing of inlet 
pressure improved the consumption of reactants and 
more homogeneous distribution. The effect of channel 
design also changed the consumption of reactants and 
consequently increases water production by Zeroual et 
al. [8], Hydration of membrane in the PEMFC is to 
maintain the performance. If it is not maintained 
properly, flooding and dehydration would occur 
which will affect the fuel cell perfonnance [9, 10]. So 
identifying the proper channel and flow field design is 
a very important task while designing the fuel cell 
which also affects the performance of fuel cell 
significantly [11]. It is clearly indicated that 
immediate attention is required for optimizing the 
simultaneous influence of operating and design 
parameters for the performance of the PEM fuel cell. 
Hence this paper has a detailed study about the 
optimization of operating pressure, temperature, 
stoichiometric ratio of inlet reactant mass flow rate 
and various rib to channel width (L:C)-1:1,1:2, 2:1& 
2:2 on serpentine flow channel of 36 cm2 active area 
of PEM fuel cell are to be studied and influence their 
performance are compared. 

II. MODEL DEVELOPMENT 

Three dimensional (3-D) PEM fuel cell model with 
serpentine flow channel of various rib to channel 
width configurations were created by Creo Parametric 
2.0 as shown in Fig. 1 



.0 (a) R:C - 1:1 

f . 

\° (b) R:C - 1:2 
A (c) R:C - 2:1 


, (d) R:C - 2:2 

Fig. 1 .Various landing to channel width ratio (a) 1:1 
(b) 1:2 (c) 2:1 and (d) 2:2 of serpentine flow channel 
of 36 cm 2 active area of PEMFC 

The modeling was done by creating all individual 
parts of the PEMFC and the dimensions of individual 
parts such as the anode and cathode GDL, solid 
polymer electrolyte membrane, the anode and cathode 
catalyst layers as shown in the Table 1. The various 
geometrical models (L: C-1:1, 1:2, 2:1 and 2:2) of 
serpentine flow channel were meshed by using ICEM 
14.5. 

After geometry modeling, the next step was 
discretization of PEMFC used by ANSYS 14.5 ICEM 
software. The Cartesian grid meshing method was 
used, which is used in the formation of hexahedral 
mesh to attain accurate results. Split block method 
used for blocking. Body fitted mesh was used and 
projection factor was set to 1. The projection factor 
determines how closely the edges of the mesh match 
up with the grid. 
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s. 

No 

Part Name 

Width 

(mm) 

Length 

(mm) 

Thickness 

(mm) 

Zone 

Type 

1 

Anode & Cathode 
Flow channel 

60 

60 

10 

Solid 

2 

Anode & Cathode 
catalyst 

0.08 

Fluid 

3 

Membrane 

0.127 

Fluid 

4 

GDL anode & 
cathode 

0.3 

Fluid 


Table l.Dimensions and Zone type, assigning of 
fuel cell 

The simulation of PEMFC was solved by 
simultaneous equations like conservation of mass, 
momentum, energy, species concentration, butler- 
Volmer equation, Joule heating reaction and the 
Nemst equation to obtain reaction kinetics. The model 
used to consider the system as 3-D, steady state and 
inlet gases as ideal condition, system as an isothermal 
and flow as laminar, fluid as incompressible, thermo 
physical properties as constant and the porous GDL, 
two catalyst layers and the membrane as an isotropic. 

A control volume approach based on commercial 
solver FLUENT 14.5 was used to solve the various 
governing equations. Three-dimensional, double 
precision and serial processing were used for this 
model. The species concentration on anode side of FL, 
O 2 , and H 2 0 were 0.8, 0, and 0.2 respectively. 
Similarly, on the cathode side were 0, 0.2 and 0.1 
respectively. The porosity at anode and cathode side 
was 0.5. Open circuit voltage was set at 0.95 V on the 
cathode and the anode was grounded. The cathode 
voltage has been varied from 0.05 V to 0.95 V used 
for solving kinetics reaction in order to get the current 
flux density, H 2 , 0 2 , and FLO fractions along with the 
flow field design. Multigrid settings were modified as 
F-Cycle for all the equations and entered termination 
restriction value was set as 0.001 for H 2 , 0 2 , H 2 0 and 
water saturation. The electric and proton potential 
values were set at 0.0001. Stabilization method 
BCGSTAB was selected for H 2 , 0 2 , H 2 0, water 
saturation, electric and proton potential. The Anode 
and Cathode reference current density was set to be 
10000 A/cm 2 and 20 A/cm 2 respectively 0.1 kmol/m 3 


was set to anode and cathode reference concentration, 
Anode and cathode exchange coefficient was set to be 
2. The Reference diffusivity of H 2 , 0 2 and H 2 0 was 
set to as 3E-5. 

Taguchi method has been used to find out the most 
optimum combination among the input parameters 
which would result in getting the maximum possible 
output which cause the performance enhancement of 
PEM fuel cell. In Taguchi method LI6 standard 
orthogonal array with 4-level and 4 factors was used 
and the parameters were considered as low, high and 
medium range values. When this orthogonal array 
was used, significance of factors and optimum 
combination can be found in 16 runs itself. The 
factors considered for the analysis were landing to 
channel ratios on serpentine flow field design (L: C- 
1:1, 1:2, 2:1 and 2:2), pressure (1, 1.5, 2 and 2.5 bar), 
temperature (313, 323, 333 and 343 K), anode and 
cathode reactants as stoichiometric ratios (S/F) of 3, 
3.5, 4 and 4.5. The theoretical value of hydrogen in 
the anode side was 4.33E-07 kg/s and oxygen in the 
cathode side was 3.33E-06 kg/s. 

III. RESULTS AND DISCUSSION 

As per LI6 orthogonal array, the inputs were given to 
the Ansys CFD Fluent analysis software and having 
all other parameters constant. The power densities for 
all 16 runs, obtained from analysis software and the 
corresponding Signal/Noise (S/N) ratios were found 
from MINITAB 17 software as shown in the Table 3. 
The rib to channel width ratio of 1:1 for serpentine 
flow field has shown maximum power densities of 
0.263 W/cm 2 and minimum power densities of 0.19 
W/cm“respectively. Similarly for L:C of 1:2 and 2:1 
having maximum power density of 0.276 W/cm 2 and 
0.199 W/cm" respectively. The minimum power 
densities for the same L:C ratios have 0.221 W/cm 2 
and 0.152 W/cm" respectively. 

For the rib to channel width ratio of 2:2 has shown 
maximum power density of 0.215 W/cm and 
minimum power density of 0.107 W/cm 2 . The 
optimization was performed for “Larger the Better” 
type of Taguchi method since power output of PEM 
fuel cell must be maximized. The S/N ratio plot for 
the same were obtained using MINITAB 17 software 
and the corresponding maximum S/N ratio gives 
better performance as analyzed based on larger the 
better as shown in the Fig.3 
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It was concluded that the design parameter such as, rib to channel ratio of serpentine flow channel having -1:2 
as L2, and the operating parameters like pressure - 2.5 bar as M4, temperature - 333 K as N2, Stoichiometric 
ratio of inlet mass flow rate - 4.5 as 04 were the optimum parameters to show the better PEMFC performance. 
The optimization results of various parameters were based on S/N ratios and the significance of each factor by 
ranking them according to their performance. Delta value of each factor available on the MINITAB 17 software 
itself was shown in Table 4. 


Run 

L:C 

Pressure 

Temperature 

Stoi.Ratio 

Power Density 
(W/cm 2 ) 

S/N Ratio 

1 

lxl 

1 

323 

3 

0.190 

-14.4235 

2 


1.5 

333 

3.5 

0.213 

-13.4159 

3 


2 

343 

4 

0.238 

-12.4805 

4 


2.5 

353 

4.5 

0.263 

-11.6141 

5 

1x2 

1 

333 

4 

0.221 

-13.0974 

6 


1.5 

323 

4.5 

0.236 

-12.5599 

7 


2 

353 

3 

0.273 

-11.2913 

8 


2.5 

343 

3.5 

0.276 

-11.1769 

9 

2x1 

1 

343 

4.5 

0.179 

-14.9336 

10 


1.5 

353 

4 

0.199 

-14.0028 

11 


2 

323 

3.5 

0.152 

-16.3448 

12 


2.5 

333 

3 

0.170 

-15.3893 

13 

2x2 

1 

353 

3.5 

0.107 

-19.4454 

14 


1.5 

343 

3 

0.128 

-17.8464 

15 


2 

333 

4.5 

0.215 

-13.3547 

16 


2.5 

323 

4 

0.196 

-14.1664 

Average S/N Ratio 

-14.096 


Table 3. Factors, levels, power density and S/N ratio for 16 runs of optimization 


Main Effects Plot for SN ratios 

Data Means 



Fig .3. Mean S/N ratio plot for L:C (L1-L4),Pressure (M1-M4),Temperature (N1-N4), Stoi.Ratio (01-04) 

The factor with highest delta value indicates higher significance factor. It was found that pressure was the 
predominant factor affecting the performance of PEMFC. The other parameters were also influencing the 
performance to a considerable extent such as, landing to channel width ratio of serpentine flow channel, 
stoichiometric ratio of inlet mass flow rate and operating temperature respectively. 
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Factors 

Level 1 

Level 2 

Level 3 

Level 4 

Delta 

Rank 

Landing to Channel 
width (L:C) 

-12.98 

-12.03 

-15.17 

-16.2 

4.17 

1 

Pressure (bar) 

-15.48 

-14.46 

-13.37 

-13.09 

2.39 

2 

Temperature (K) 

-14.37 

-13.81 

-14.11 

-14.09 

0.56 

4 

Stoichiometric ratio 

-14.74 

-15.1 

-13.44 

-13.12 

1.98 

3 


Table 4. Mean S/N ratios, Delta and Ra nk for each level of factors 


The percentage contribution of individual parameters, P-test and F-test on the serpentine flow fields for the 
perfonnance of PEMFC has been shown in the Table 5. It has been observed from the Table 5, operating 
pressure has been contributed to be 43.5 % , operating temperature was 0.72 %, the stoichiometric ratio of the 
reactants and L:C has contributed 8.72 % and 26.3 % respectively of the PEMFC perfonnance 


Factors 

DOF 

Sum of squares 

Variance 

F-test 

P-Test 

Contribution (%) 

Pressure 

2 

0.006704 

0.003352 

88.71 

0.283 

43.50 

Temperature 

2 

0.000043 

0.0000215 

27.26 

0.566 

0.72 

Stoichiometric ratio 

2 

0.001466 

0.000733 

19.39 

0.825 

8.72 

L:C 

3 

0.006171 

0.002057 

5.4 

0.036 

26.30 

Pressure & 
Temperature 

1 

0.000074 

0.000074 

0.22 

0.428 

0.02 

Pressure & L:C 

3 

0.003672 

0.001224 

0.81 

0.059 

15.24 

Error 

2 

0.000151 

0.0000755 



5.50 

Total 

15 

0.035937 

0.007537 

141.8 

2.197 

100 


Table 5.The percentage contribution of individual parameters of serpentine flow channel 


Also the combined effect of combination of pressure 
with temperature and pressure with L:C has shown 
0.02 % and 15.24 % respectively contributing to peak 
power performance of the PEMFC. 

CONCLUSION 

The combined effect of all the parameters exhibited a 
different response compared to their individual 
effects. The maximum power density of optimizing 
the four different parameters on serpentine flow 
channel of 36 cm2 active area of PEM fuel cell using 
Minitab 17 provides 0.276 W/cm 2 from L:C-1:2 with 

2.2 bar operating pressure, 343 K temperature and 3.5 

2 

stoichiometric ratio of inlet reactant gases and R 
value was arrived 96.9 %. The effect of operating and 
design parameters was affecting the performance of 
PEM fuel cell considerably. 


REFERENCES 

[1] Manso, A. P., Garikano, X., Mujika, G. M., 
Marzob, F.F., Barranco, J. "Influence of geometric 
parameters of the flow fields on the perfonnance of a 
PEM fuel cell. A review."International Journal of 
Hydrogen Energy. 37, pp. 15256-15287. 2012. 

[2] Shimpalee, S., Greenway, S., Van Zee, J. W. "The 
impact of channel path length on PEMFC flow-field 
design." Journal of Power Sources. 160, pp. 398-406. 
2006. 

[3] Lakshminarayanan V, Karthikeyan P, 
Muthukumar M, Senthilkumar A P, Kavin B, 
Kavyaraj A, ‘Numerical investigation of performance 
studies on single pass PEM fuel cell with various flow 
channel design’, Applied Mechanics and Materials, 
Vols. 592-594 , pp 1672-1676, 2014. 

[4] Siefert, N. S., Litster, S. "Voltage loss and 
fluctuation in proton exchange membrane fuel cells: 


485 

IJTSRD | May-Jun 2017 
Available Online @www.ijtsrd.com 
































International Journal of Trend in Scientific Research and Development, Volume 1(4), ISSN: 2456-6470 www.ijtsrd.com 


The role of cathode channel plurality and air 
stoichiometric ratio." Journal of Power Sources. 
196(4), pp. 1948-1954.2011. 

[5] Lakshminarayanan V &Karthikeyan P. 
Optimization of Flow Channel Design and Operating 
Parameters on Proton Exchange Membrane Fuel Cell 
Using Mat lab. PeriodicaPolytechnica Chemical 
Engineering.Budapest Univ Technology Economics. 
2016, 60, 3; 173-180. 

[6] Khazaee, I, Ghazikhani, M &Mohammadiun, M , 
'Experimental and thermodynamic investigation of a 
triangular channel geometry PEM fuel cell at different 
operating conditions', Scientialranica, vol. 19, no. 3, 
pp. 585-593,2012. 

[7] Yan WM, Chen CY, Mei SC, Soong CY & 
Chen F. 'Effects of operating conditions on cell 
performance of PEM fuel cells with conventional or 
interdigitated flow field, Journal of Power 
Sources,2006, 162, no. 2; 1157-1164. 


[8] Zeroual,M.; BelkacemBouzida,S. ; 
Benmoussa,H.;Bouguettaiaa,H.(2012). Numerical 
study of the effect of the inlet pressure and the height 
of gas channel on the distribution and consumption of 
reagents in a fuel cell (PEMFC), Energy procedia, 18, 
205-214. 

[9] Xianguo Li, Imran Sabir&Jaewan Park. (2007) „A 
flow channel design procedure for PEM fuel cells 
with effective water removal”, Journal of Power 
Sources, 163, pp 933-942. 

[10] Chi-Young Jung, Chi-Seung Lee & Sung-Chul 
Yi.(2008) Computational analysis of transport 
phenomena in proton exchange membrane for 
polymer electrolyte fuel cells”,Journal of Membrane 
Science,309,ppl-6. 

[11] AtillaBiyikoglu. Review of proton exchange 
membrane fuel cell models, International Journal of 
Hydrogen Energy,2005, 30, 1181 - 1212. 


IJTSRD | May-Jun 2017 
Available Online @www.ijtsrd.com 


486 


